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Abstract 

Using gradient expansion method the particle number, the energy and the action 
were calculated at the formation of a vortex-skyrmion. Contrary to other papers on this 
subject no approximation of single Landau level was used. Taking into account the 
nearest Landau levels essentially changes the expression for the vortex energy and gives 
some simple physical interpretation of the main results. It is shown that the formation 
of a vortex skyrmion gives the gain in the thermodynamic energy and therefore they 
must be spontaneously created near the odd fillings. The Hopf term in the skyrmion 
action is calculated and corresponds to a fermionic picture. 

A large number of theoretical works is devoted to the consideration of macroscopical 
skyrmioniclike excitations near odd Landau level fillings in 2d systems under Quantum Hall 
Effect conditions. In the work [1] the phenomenological approach of the Chern-Simons theory 
was used. It was shown the existence of such exitations and their energy was calculated. 
In [2] the energy was calculated using Hartree-Fock approximation for single Landau level 
projected wave functions. The gradient expansion method was used in [3] for the same model 
assuming the large size of a skyrmion. The energy and the particle number were found 
analitically. These results were confirmed in [4] and the general technique was developed 
for the calculation of arbitrary terms in gradient expansion. In the model of projected 
functions the skyrmion action was calculated [5] with the topological term (see the proper 
discussion [6]). The experiments give some evidence of skyrmion formation [7] but some 
experimentalists cast some doubts concerning these results[8]. 

Single Landau level projected functions approximation is usually justified by the large 
value of the cyclotron energy Tiuo c compare to Coulomb interaction e 2 /ln where the magnetic 
length l 2 H = ch/eH and H is the external magnetic field. The appropriate calculations are 
rather elaborate and contain some additional assumptions. The final expressions for the 
electron energy and the density are obtained by cumbrsome calculations without any direct 
physical interpretation. In the present work we show that the approximation of the projected 
wave functions is not sufficient for the adequate skyrmion description. Some preliminary 
results were published in [9] and the present paper is more complete with the justification 
of the results by gradient expansion method. 

Skyrmions correspond to the rotation of the second quantized operators-spinors using 
nonuniform rotation matrix U(r). In that case the initial spinors in the laborator frame i/j 
are expressed in terms of the new spinors x with spin up in the local frame according to the 
relation if)(r) = [7(f)x(f).The matrix U(r) can be parametrized by three Euler angles 
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U z {a) = cos — + 1 sin —<r z 

Uy(P) = cos - + i sin -o- y , 

where a x , a y , o~ z are Pauli matrices. At large distances from the core for finite value of g-factor 
the average spin must be oriented along the external magnetic field to have a finite value 
of the total energy. Therefore the angle (3, counted from the magnetic field direction must 



rapidly vanish for r — > oo (it can be shown exponentially) We assume that the matrix U(r) 
has no singularities at any r.That corresponds to a nonsingular behaviour of the matrices 



dxi 



where k = x,y and I = x,y,z, Pauli matrices 
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The expression for the functions Q l k are easily obtained by the direct differentiation of U 
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The non trivial topology defined by the matrix U(r), is connected with the properties of the 
mappings a (r), 7(f), where r is varied along a circumference of large radii. The degree of 
the mapping of 2d plain x, y on the sphere parametrized by the angles 7 and (3 coincides 
with the degree of the mapping of the circumference on the circumference i.e. the winding 
number for the vortex singularity of 7(f). For a nonsingular f^(r) it is necessary to have 
coinciding singularities in 7(f) and a(r) at the the point where cos/5 = —1. Therefore in 
the matrix U must be represented all three Euler angles and the corresponding spinor %p{r) 
has vortex singularity of it's phase at large distances with the integer winding number and 
is simple in 2d plain. For that reason it is more appropriate to talk about a non singular 
vortex with skyrmionic core by the anology with H e^. The winding numbers are here any 
integgers contrary to the case of He^ where they are even [10]. The integral 
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rottt z d 2 r = Q 



is a topological invariant and is expressed directly through the change of the phase for the 
spinor ip along a large circumference. 

The pair interaction Lagrangian for electrons : 
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dip 1 
dt 2m 



d 2 rdt + 



V(r — r')i[) + (f)^ (r')il)(r')i/;(f)d r'd rdt — gH l ip + o z ipd rdt 



(2) 



transforms by the matrix U into the Lagrangian in terms of spinors \ 
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d rdt 



- - J V{r - r')x + (r)x + (^)x(f')x(r}d 2 r'd 2 rdt + gH J U + a z U)xd 2 rdt (3) 

without any approximation. We assume that matrix U depends on the time also. The 
quantity Vt\ corresponds to (jl]) but with the time derivative. A Q is the vector potential of 
the external uniform magnetic field. The Lagrangian (§) can be easily obtained by direct 
differentiation taking into account some identity obtained by the differentiation of the of the 



The size of the vortex-skyrmion core is defined by Coulomb interaction which tend to 
increase the domain with large change of the density and by Zeeman energy which tend to 
diminish this domain due to the positive energy of the spins non aligned to the external 
magnetic field. At large distances Coulomb is small and only Zeeman energy and the energy 
connected with spin direction gradients are essential. That gives the exponential decay of 
non aligned spins. For a small value of the g-factor the size of the core will be large and the 
derivatives of the matrix U become small that justify the use of the gradient expansion for 
the calculations. We shall treat this case interesting mostly in topologically invariant terms 
not dependig on deformations of U. That allow us to concider no Zeeman energy or taking 
it into account only in the first order of the perturbation theory. 

The effective action depending on the matrix U is defined by the integration out the 
fermions in the correspondent expression for the statistical sum and has the form [11] S — 
iTr In G, where G is electron Green function. The Trace is over whole set of variables 
including spin and space-time. This action will be calculated by the gradient expansion in U 
using Hartree-Fock approximation which we justify later. 



1 Green function in Hartree-Fock approximation 

We use Hartree-Fock approximation for the calculation of Green function which is valid for 
fully filled Landau level with the accuracy of the order Vi nt /hu c assumed to be small. Such 
an approach allows to write the electronic Hamiltonian as 
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The interaction term becomes the sum of exchange and direct interaction and small Zeeman 
energy is ommited. In the direct unteraction we must ommit zero Fourrier component for 
Coulomb case due to compensating uniform charge (see e.g. [11]). 

We assume that all the average quantities are close to their uniform values for fully filled 
Landau level. We restrict our consideration by the model of the local exchange neglecting the 
difference r — r 1 that formally corresponds to the interaction distance in V(r) small compared 
to the magnetic length. Probably the transition to the local exchange is quite unessential 
and topological terms does not depend on it but the use of such an approximation greatly 
simplifies the calculations. The final results are easily generalized on the case of the nonlocal 
exchange. Therefore we shall deal with the hamiltonian 
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d 2 r (5) 



where 7 is the exchange constant and Vo is the constant of direct interaction wich we shall 
omit as in the case of Coulomb interaction, p is the average density. We suppose that 
X spinors are oriented in z-direction coinciding with the average local spin direction in the 
laborator frame. Assuming a small value of the angle velocities Q l and their derivatives dQ ~ 
Q 2 we shall calculate Green function by the perturbation theory separating the hamiltonian 
H = H + Hi + H 2 into three parts 
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H 2 = — I \ ' 
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dr k 



X d 2 r (8) 



We use grand canonical ensemble introducing chemical potential \i. 

Green function for the unperturbed Hamiltonian is Green function for noninteracting 
electrons in uniform magnetic field. We assume that at the last filled Landau level only the 
lowest spin sublevel is occupied. For simplicity we consider here only the case when Landau 
index s = when unperturbed Green function acquires the form 

G (fy,t,t') = -i < T X (r,t) X + (f*,t') >= W ^( W ) e ^'-')$ sp (f)$: p (0^ (9) 

p,s J 27r 

Here T is a symbol of time product for fermions, the sum is over all s and p, spin matrices 
g s {w) corresponds to fully occupied Landau level s = with spin up in the local frame. Other 
states are empty. All quantities are measured in units where Hq = 1, Ijj = 1 and h = 1. The 
normalized wave functions $ sp (r) corresponds to eigenfunctions in Landau gauge. 
It is easy to find matrices g s (u) for the hamiltonian (|j) 

0o(w) = —— ( * - (10) 

9s(uj) = — s (11) 

where 5 — > (+0), the energy of zero Landau level l/2m is included in chemical potential. 

Using the expression (^) we can develop the perturbation theory expansion for Green 
function of the total hamiltonian G = Gq + G% + G 2 + .... The corresponding diagrams are 
shown in fig 1 up to the third order. 

It is usefull for further calculations to represent the operator entering (|7]) in the form 

-n l k(Tl (-id k - A ok ) = -(Q l + 7T- + n l + n + )a t (12) 
m m 



The operator 7r + $ sp = y2(s + l)$( s+ i) p increases Landau index,and the operator ir $ sp = 

V / 2s$( s -i)p decreases Landau index what follows directly from the properties of oscillator 
wave functions and 

n l x y Q l = _x y 13 

+ 22 y J 

The correction of the first order in Q to Green function is 

G x (fy,t,lf) = J e-^-^g^^^in) 

x [I(^(f 1 ) 7r - + QL(n)ir + ) + n\{r x )](Tig s >{u') 
m 

x e^^^^^^ndt, (14) 

We are interested in Green function G(r,r;t,t + 5) defining the density. In that case the 
only nonvanishing terms are with s = 0, s' = 1 , s — 1, s' = and s = s' = 0. The others 
vanish because of analytical properties of g s for s > or have some additional derivatives. 
It is possible to neglect the difference t± — t for the same reason. We do not consider at first 
the term with fl\ which can be easily calculated. Integrand in flT4]) is vanishing rapidly when 
\r — fi\ is increased above magnetic length. Therefore it is possible to expand fi(ri) in terms 
of r — T\ up to linear terms. Performing the summation in p we get the following expression 
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Here we introduce the functions 



2tt 



Rooi^n) = — exp( 

Z7T 



<9r fc 



■d 2 n 



^i(r,fi) = i-[l-£-^ 

Z7T Z 



Rio(n,f) 

Roi(ri,r) 



1 xi - x - - y) 

1 x - xi - i(yi - y) 
2tt 



exp 
exp 
exp 



(r 




2 




4 




(r 


-fi 


) 2 




4 




(r 


-n] 


2 




4 




(r 




2 


4 



(15) 

(16) 
(17) 
(18) 

(19) 



After the integration in r*i the term wiothout the derivatives of Q vanish because the inte- 
grand is an odd function. Doing all integrations in([L^) we get using (|13|), (0), (|10l ) 
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here we consider only the main spin independent contribution to g\ ~ m. In the term with 
Vl\ we must consider only the case s = s' = others are of the order of m and can be 
neglected. Only the terms with the first order poles in u are essential that results in 
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Finally the first order correction to Green function acquires the form 

Gi(f, f, t, t + 5) = ^- \(cr z (Ti - aia z )divtt l - z'((l + cr z )a l + (7,(1 + a z ))rotQ l 

-(divVt - irotVt ) - iQr t - 
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For the density p = —iSpG(r, r,t,t + 5) we obtain 

p(r,t) = -'-(l -rotQ z 
2n 
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This result was firstly obtained in the approximation of the projected wave functions [3]. 
In the terms of unprojected wave functions this result has a simple physical meaning: the 
lowest Landau level in the local effective field H e ff = 1 — rotQ z is fully filled and the density 
coincides with the local density of states p = l/2nl 2 (H e ff). 

The results for the first order term in G contain the derivative of Q that corresponds to 
second order in Q. Therefore it is necessary to calculate the formally second order term in G 
which can be written down simbolically as 



G 2 (f, r,t,t + S) = ^rG (tt l + 7r- + ^tt+V^o^+tt- + tfLn + )<T V G + G H 2 G 
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In the second order we can neglect the additional derivatives of Q. Because the last term 
is formally of the second order we must take in it the states with s = with the second 



contributions with gogigo and gigo9i, where g corresponds to one of the three Green functions 
in this expression. The integration over space variables and p are easily performed and we 
get 

G 2 = — — I {goaig^gott 1 ^ + g 1 a l g a l ,g 1 ttin l +]e i " s ^d 2 rdt 

m 271 J 2lX 

To calculate the second order correction to the density we need trace of this expression. 
Transposing the matrices under the trace we get in the first case the square of diagonal 
matrix g^uS) having the second order pole in u. Using the standart expression for the 
proper integral we obtain 



SpG 2 = J [-2glQ} + Q}_ + 2c^fi , _ 



d 2 rdt = 



Therefore the density is given by the expression fl2T|) up to the second order. 

2 The calculation of the action, the energy and the par- 
ticle number 

For the calculation of the action S = iSplnG we use the perturbation theory for Green 
function and expand the logarithm 

S = iSp[H 1 G + -HxG Q HxG + H 2 G + -(H 2 G HiG + HiG H 2 G ) + -^HxGqH 1 GqH 1 Gq+ ..] 

Here we omit nonessential part without the rotation matrix and only terms up to the third 
order are retained. In this section we limit our consideration by the second order terms. The 
action of the first order is given by 

S x = iSp I n\a igo (uj)e^ 5 ^p + iSp f a m {uo)e^{n l + 7x~ + n l _ir+)% p {r)% p {r)d 2 rdt 
J Ztt 2tt J 2tt 

The second term can be calculated using integration by parts and we get 

S 1 = -—J Q z t drdt - — — J (idivVt z - rotQ z )d 2 rdt (22) 
The action of the second order consists of two terms one of which contains two Hi H^. 
S' 2 = ^Sp J (n l _7r + + n l + n-)a l9s (uj)$ sp (f>)$: p (f>) 

(^7r + + Q^a.gA^s'pin^v^e^d'r'd'rdtp 
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In this expression we make no difference between Q(f,t) and f2(r 1 ',t / ) neglecting derivatives 
of fi.The only nonvanishing terms correspond to s = 0, s' = 1 and s = 1, s' = 0. Using the 
properties of the operators tt + ,tt~ and performing simple integration it is easy to get the 
second order action including the second term with H 2 
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The symmetric part with I — I' of the first term is compensated by the proper part of the 
second term. The antysimmetric part containes only 1,1' = x,y and can be rewritten using 
the identity rotfl z = 2(Q^.Q^ — VtyVt y x ) which is obtained by differentiation using the identity 
U + U = 1. Finally we get the second order action 
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The total action up to the second order inclyding ( p2] ) is 
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Here we add the term with chemical potential. 

These results are in full agreement with the authors preliminary publication [9] and have 
simple physical interpretation. The electrons locally fill the the lowest spin sub level in the 
local "effective" magnetic field H e ff = 1 — rotQ z . They fill it fully so the electron density 
coincides with the local density of states. This circumstanse justifies the use of Hartree-Fock 
approximation in the same manner as it does for the fully filled Landau level. The "effective" 
magnetic field is less than the external one for vortices with positive Q. It must be stressed 
that the gap beween spin sublevels is defined by the exchange. The electrons with the spin 
direction opposite to average spin direction will see the "effective" field H e ff = 1 + rotfl z . 

Let us consider more detailed expression for the energy in a stationary case taking into 
account the change in the interaction terms due to the variation of the density. For the 
Coulomb interaction one must assume the existence of the compensating positive background 
recquiring to vanish the zero Fourrier of the interaction. Because of that there is no linear 
term in the density variation for direct interaction. It's expansion begins from the second 
order 

1 f e 2 .x„,^,d 2 rd 2 r' 
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The exchange energy begins with a linear term in the variation of the density. For it's 
calculation one must express the density of the exchange energy as a function of the magnetic 
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field for fully filled Landau level — differenciate it in magnetic field and multiply by 

the variation in effective magnetic field due to vortex formation 

r 3e 2 
E ex = / — ^^rottt z d 2 r 
J AV2l^l H 

Also it is necessary to take into account an energy correction due to nonuniform direction 
of the average spin direction J J(^-) 2 d 2 r where the quantity J = —k== j-, (see e.g. [3,4]. 
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and add Zeeman energy. Finally we get for the change of the thermodynamic energy due to 
vortex-skyrmion formation 

e 2 f rotQ z (f)rotQ z (r f ) >9 >9 . [A ^,dni.~ 1 , , 9 . „. 
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We assume here the value of chemical potential fi = %r corresponding to the fully filled 
lowest Landau spinsublevel at large distances from the core. By assumption the cyclotron 
is large compare to the interaction energy and consequently the quantity F is negative for 
positive topological number Q defining the change of the total number of electrons SN = —Q 
at vortex formation. Therefore the formation of nonsingular vortices-skyrmions give the gain 
in the thermodynamic energy and they must be spontaneously created. 

It is possible to calculate also the energy of one particle exitations. This can be done 
by the investigation of the poles for G(r,r,u). However it is more simple to use the fact 
that these energies corresponds to to the variational derivatives of the total energy over the 
local density of the electrons with the different spin direction. It is possible to write down 
the quantity F as a functional depending on the densities p+,p~ where signs denote the 
direction along and opposite to the average local spin 
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In this expression 7' ' p + = y f — ^f^y f ro ^ Z an d </'p+ = ^- Variation in densities 
p+,P- gives us the energy of the hole in the filled states 

5F 2e 2 /tt, 3. 



£ » = -^ = 77V§( 1 -i'> t ^ + (27 > 

and the energy of the electron with the locally reversed spin 

All previous results are obtained by using the perturbation theory in f2. That means that 
one particle energy corrections due to f2 must be small compared to the initial gap separating 
the lowest spin sublevel from the higher one. Therefore the largest additional term in one 
particle energy %uo c /2l 2 H rotVL z ~ %uo c /2l 2 H / L 2 where L c is the core size must be small compare 
to the exchange energy ~ e 2 /l H . That gives the inequality L c 3> lH\Jhu c /e 2 . The core size 
is defined by the competition of Coulomb energy e 2 /L c and Zeeman energy gHpL 2 that 
gives L\ ~ (e 2 l 2 H )/ gH ^> P h [{TujJ c Ih) I 'e 2 } 3 ^ 2 '. Therefore for the validity of the obtained results 
g-factor must be small enough. The addition of the electron with reversed spin increase the 
thermodynamic energy due to exchange and decrease it by Coulomb interaction for positive 
Q other terms are essentially less due to small value of Q. Coulomb energy increases with Q 
and inverse size of core whereas exchange energy is local and does not depend on the core 
structure. Therefore it is possible to have metastable bound states and even the states below 
chemical potential. In that case the total amount of the electrons with "wrong" orientation 
in the external magnetic field will be diminished by one compare to the vortex without bound 
electron. The thermodynamic gain by the creation of the vortex will be conserved. For real 
experimental conditions the cyclotron energy is not very large compare to Coulomb energy 
and g-factor is also not very small and therefore these statements need computer check. 



3 Hopf invariant in the action 

Besides the mapping degree Q there is also topological Hopf invariant corresponding to the 
linking number of the lines with constant values of n(r, t) for nonstationary matrix U (see 
e.g. [13]). This invariant must be contained in the action with the coefficient determining 
the statistics of vortices or more precise the phase change when they are transposed [13]. 
This coefficient was calculated in [5] in the model of projected wave functions. However this 
calculation contains some additional assumptions and was questioned in the discussion [6]. 
The authors of [6] reference to their previously made semiclassical calculation which use the 
frequency and the momentum as a good quantum numbers for Green function calculation 
and cite the proper result. In this section we make a calculation of this coefficient directly 
in the limit of large magnetic field without additional assumptions. For that one needs to 
find the action in the third order of Q including only the terms containing two space and 
one time component of Q l . 

Third order perturbation terms in the action are 

S = iTr^H.G^Gv + (H^o^Gq) + ^H^H^H^] (29) 

The terms up to the second order were calculated in the previous section. However time 
derivatives were not included in these terms what is essential in the third order calculations. 



We begin with the second order in H\ term containing VL l t omited in the previous section 
assuming Q t ~ After an easy calculation we get 
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With the needed accuracy there are one term with s = s' = and two terms with s — 1, s' — 
0; s = 0, s' = 1. Introducing the new integration variables R = (r + r 1 )/2,p = r [ — f and 
expanding fif(r), f2(r*i) in p up to the first order we obtain for the first term 
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Using expressions (|T0|) , (]TT|) and performing the proper integrals in u and p we have 
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The terms with s — 0, s' — 1; s = 1, s' = are calculated in the same way 
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There is also the term with two Q l considered in the previous section not taking into 
account time derivatives. We can proceed quite similar introducing the new integration 
variables T = (t + t')/2,r = if — t and expanding in r. With the required accuracy it is 
sufficient to consider omly the terms with s — 1, s' — and s — 0, s' — 1 others contain 
extra derivatives. Omiting the second order term without time derivatives we get 
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Representig r as the derivative of the proper exponent in the integrand and integrating by 
parts we have 
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In order to compare various terms in action they must be expressed in some standard 
way through Q l and it's derivatives. Using the identity d t VL l k — d k VL l t = 2e^ m Vt , t VL k n where e Zjm 
is antisymmetrical tensor it is easy to transform the previous expression to the following 
form 

Sf = -\j (pf x VQ l t + e ljrn Q l x Q m Qi) (33) 
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Among the third order terms the simplest is containing H2 and Hi 
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The other terms contain extra derivative. Only the term with the divergence gives nonvan- 
ishing contribution 

The last remaining third order terms contain only H\ in the combination with two space 
and one time fl l 
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Only the terms with s = s 2 = 0; Si = 1 and s = s 2 = 1; s\ = gives the needed third order 
and no derivatives of Q must be taken into account. 

The calculations in the first case are very similar to the calculation of fl3T|) , (|35| ) for double 
pole in the production of go- Using the standard form in terms of Q\ and Qi we get 



i^a Jl01a / + \ E / n?(^') 2 ^ (36) 

The calculation of the case s = s 2 = 1; s% = is standard and gives 

sl = -\J ^f d ^+W J mnr-^n^+ypa^ J 

(37) 

The final expression for the part of the action containing the Hopf invariant is obtained 
by the summation of all found terms fl3?D , (|36|), fl35|) , (^3|), fl32"D, Some points in this 

summation must be explained. The terms with the factor I/7 in S3, S 2 , S 2 i combine to the 
total derivative 

4 ^ E / HdMnJtf) - v x (n l t n 1 )}^ = o 



because fr for I 7^ 2 tends exponentially to zero at large distances. In the same way the 
term containing the other total derivative 

d 2 rdt 



/— * — » ci v fit 

[n l t divn m + (n m v)n l t }—- = o 
Z7T 
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vanish. 

The terms with Qf(£l 1 ) 2 and Q l t (Q z Q l ) combine also to zero after some algebraic manip- 
ulations with S 2 . Finally only antisymmetric terms in S3 are retained 

d2rdt , 1 /"niw „ Xl d2rdt f JlmrJrtl „ ^m^^ 



2ty 



s 3 = Ismwh i n l t n h x + i /^v x n'— - / e^n' x n 

4 j 2n 2 J 2n J 

Calculating the trace and using the identity for V x fr we get 

S 3 = e ljm [ n\W x (38) 

J Z7T 

According to [13], [12] the Hopf invariant can be expressed in terms of Q 

H=-L [ e ljm nln j x n m d 2 rdt 



Therefore the Hopfian term in vortex action is 

S H = tiH (39) 

. This result coincide with those in [6] . According to the accepted point of view that means 
[13] the fermionic nature of vortices-skyrmions. 

We must remain that the Hopf invariant H is integer for the mapping of the sphere S3 
on the sphere S2 that requires in our case vanishing of Q at large f, t i.e. Q = 0. If the 
vortices are present at any time and Q 7^ then the topological invariant corresponding to 
the linking number changes a bit and is contained in some part of H [15], [16]. However the 
found relation between S and H holds on as well as the statement about fermionic nature 
of vortices. 

4 Conclusions 

In this work the theory of nonsingular vortices-skyrmions was considered for 2d electron 
systems in a strong magnetic field. It was assumed that g-factor is sufficiently small to 
justify the perturbation theory in the derivatives of the rotation matrix. No other assump- 
tions beside the large value of cyclotron energy Tiuj c compared to Coulomb energy e 2 /lu was 
used. It was shown that currently used approximation of the single Landau level projected 
wave functions is not sufficient for the adequate description. Taking into account the near- 
est Landau levels is essential and changes the thermodynamic energy for vortex-skyrmion 
formation diminishing it's value by —hu c /2 large by assumption. This fact must lead to 
the spontaneous creation of vortices with positive degree of mapping near odd Landau level 
fillings. Simple physical picture emerges for small enough g-factor providing validity of the 
perturbation theory : inside the vortex core exists an additional effective magnetic field with 
the total flux containing the number of the flux quanta equal to the degree of the mapping. 
The lowest Landau spin sub level in the total (external plus effective) local field is fully filled 
and is separated by the exchange gap from the higher energy states. The total number of 
electrons changes because of the change of their density coinciding with the local density 
of states. As aconsquence the electrical charge of the core emerges which equals tjj that 
of electrons expelled from the core in the number of flux quanta for a positive degree of 
mapping. 

There is a possibility to have bound electron states in the positive charge of the core. The 
position of the proper energy is defined by negative Coulomb energy and positive exchange 
energy not depending on the core. Under made assumptions it is above chemical potential. 
In real experimental situation this statement must be checked numerically. 

The topological Hopf term in the vortex-skyrmion action is found. The calculated coeffi- 
cient before the Hopf invariant coincides with semiclassical result [6]. Therefore the vortices 
are fermionic-like. The properties of vortices-skyrmions are close to "composite" fermions 
introduced phenomenologically. in [14]. 
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